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We investigate the space of quantum operations, as well as the larger space of maps which are 
positive, but not completely positive. A constructive criterion for decomposability is presented. A 
certain class of unistochastic operations, determined by unitary matrices of extended dimension- 
ality, is defined and analyzed. Using the concept of the dynamical matrix and the Jamiolkowski 
isomorphism we explore the relation between the set of quantum operations (dynamics) and the set 
, of density matrices acting on an extended Hilbert space (kinematics). An analogous relation is es- 

■ tablished between the classical maps and an extended space of the discrete probability distributions. 
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(SJ . I- INTRODUCTION 

The theory of quantum information has been studied for at least forty years - consult e.g. some early papers 
by Ingarden and his group written in sixties and seventies 0, 13 HM EH- However, rapid progress of this theory 
occurred only in the last decade 0, 0, H3, H3, 0, due to a synergetic feedback with several recent experiments 
H ' on quantum physics, motivated by information encoding and processing. The possibility to encode the information 
7-H \ in quan tum states triggered an increasing interest in the structure and the properties of the set of quantum states 
O ■ EM E1E1 lli Usually one considers states acting on a finite dimensional Hilbert space, TLn- 

To characterize the way that information is processed according to the laws of quantum mechanics, one needs to 
describe the dynamics of the density matrices. In many cases it is sufficient to consider discrete dynamics, which 
maps an initial state p into the final sta te p 1 . Such maps are often called quantum channels, and their theory is a 
Qh! subject of considerable recent interest p. fl3L l45t l7q . 

(Jj ■ The set of quantum states Ai^ N ^ consists of density matrices of size N, which are normalized, hermitian, and 
positive definite. A quantum channel is called positive, if it maps the set of positive operators into itself. However, 
since a physical system under consideration may be coupled with an environment, the maps describing physical 
processes should be completely positive (CP), which means that all their extensions into higher dimensional spaces 
remain positive. Any CP map may be (not uniquely) represented in the so-called Kraus form, by a collection of 
Kraus operators |60| . 

Even though positive, but not completely positive, maps cannot be realized in the laboratory, they are of a g reat 
theoretical importance, since they can be used to detect entanglement in density matrices written on paper |37l l4l|. 
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The phenomenon of quantum entanglement seems to be crucial in the theory of quantum information. The structure 
of the set of all completely positive maps, CVn, is relatively well understood [H|, but the task to characterize the 
larger set of all positive maps is far from being completed |6l], llH |96| . 

The aim of this work is twofold. On one hand we present a concise review of recent development concerning the 
properties of the set of quantum maps. On the other hand, we present several new results in this field. In particular, in 
section 2 concerning the matrix algebra, we define a simple transformation of a matrix, called reshuffling, and establish 
a useful lemma: The Schmidt coefficients of a matrix A, treated as an element of a composite Hilbcrt-Schmidt space 
of operators, are equal to the squared singular values of the reshuffled matrix A R . In section 3 we investigate the 
properties of the dynamical matrix D associated with any linear map |94j - it may be obtained by reshuffling of the 
superoperator L which describes the map. Furthermore, we define a class of unistochastic maps, which are determined 
by a single unitary matrix U of size A^ 2 . 

For any CP map the corresponding dynamical matrix D is hermitian and positive definite. Making use of the eigen 
representation of D one may define the canonical Kraus form of the operation. In section 5 we present a comparison 
between different classes of quantum maps, acting on the space of density matrices, with the classes of classical maps 
acting on the simplex Ajv-i of discrete, A^-point probability measures. For concreteness we present in section 6 
certain exemplary maps acting in the space of qubits - the density matrices of size N = 2. 

In section 7 the positive, co-positive and decomposable maps are analyzed and a constructive criterion for decom- 
posability of a map is provided. Furthermore, for any quantum map we introduce three quantities useful to locate 
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it with respect to the boundaries of the sets of positive (CP, CcP) maps. In the subsequent sections we explore the 
Jamiolkowski isomorphism j^, which relates the set of quantum maps acting on Ai^ N \ with the states from Ai^ N \ 
which act on the extended Hilbert space 7ijv <8> Ti-N- This isomorphism, formulated as well in the quantum, as well as 
in the classical setup, allows us to discover an analogy between objects, sets and problems concerning quantum maps 
and quantum states. This key issue of the work explains its title: there exists a kind of duality between between 
properties of the set of quantum maps (dynamics) and the the set quantum states of a composite N x N system 
(kinematics). 



II. ALGEBRAIC DETOUR: MATRIX RESHAPING AND RESHUFFLING 



In this section we are going to discuss some simple algebraic transformations performed on complex matrices, 
which prove useful in description of quantum maps. In particular, we introduce a convenient notation to work in the 
composite Hilbert space Hn £3> T~Cm or in the Hilbert-Schmidt (HS) space of linear operators, Tins- 

Consider a rectangular matrix Aij, i = 1, ...,M and j = l,....N. The matrix may be reshaped, by putting its 
elements row after row in lexicographical order into a vector dk of size MN, 

d k = A i:j where k = (i-l)N + j, i = l,...,M, j = l,...N. (2.1) 

Conversely, any vector of length MN may be reshaped into a rectangular matrix. The simplest example of such a 
vectorial notation of matrices reads 



A 



An A 12 
A21 A22 



{A u , A u , A 21 , A 22 }. (2.2) 



The scalar product between any two elements of the HS space 7ins (matrices of size N) may be rewritten as an 
ordinary scalar product between two corresponding vectors of size -/V 2 , 

(A\B) = Tr A^B = d* -b=(a\b). (2.3) 

Thus the HS norm of a matrix is equal to the norm of the associated vector, = |a| 2 . 

Sometimes we will use both indices and label a component of a by ay. This vector of length MN may be 
linearly transformed into a' = Ca by a matrix C of size MN x MN. Its elements may be denoted by Ckk 1 with 
k, k' — 1, MN, but it is also convenient to use a four index notation, Cma where m, n = 1, TV while fi, v = 1, M. 

nv 

In this notation the elements of the transposed matrix are C„\ u = C nv , since the upper pair of indices determines 

nv m / i 

the row of the matrix C, while the lower pair determines its column. The matrix C may represent an operator acting 
in a composite space H — Hn <8> "Hm- The tensor product of any two bases in both subspaces provides a basis in H., 
so that 



Cmu = {e m ® f a \C\e n ® f v ), (2.4) 

nv 

where the Roman indices refer to variables of the first subsystem, Ha = Ti-N, and the Greek indices to the second, 
TLb = Hi/. The trace of a matrix reads TrC = Cmn, where summation over the repeating indices is assumed. The 

mfj, 

operation of partial trace over the second subsystem produces the matrix C A = Tr^C of size N, while tracing over 
the first subsystem leads to a M x M matrix C B = Tr^C, 

C mn — Cmfi, and = Cmfi- (2.5) 

n/t ^ rav 

If C is a tensor product, C = A® B, then Cmp, = A mn B a v 

nv 

Consider a unitary matrix U of size A^ 2 . Its iV 2 columns (rows) Uk, k = 1, iV 2 reshaped into square matrices Uk of 
size N form an orthogonal basis in Tins- Using the Hilbert-Schmidt norm, ||A||#s = (TrAAt) 1 / 2 , we normalize them 
according to Ak = Uk/\\Uk\ \hs an d obtain the orthonormal basis, (Ak \ Aj) = Tr A\Aj = Skj ■ Alternatively, in a double 
index notation with k — (m — 1)N + [j, and j = (n — 1)N + v this orthogonality relation reads (A m ^\A nv ) = 5 mn 5nv 
Note that in general the matrices Ak (also denoted by A m>1 ) are not unitary. 

Let X denote an arbitrary matrix of size A^ 2 . It may be represented as a double (quadruple) sum, 

N 2 N 2 

l X > = E E C ^\ A k) ® \A 3 ) = Crn^\A m ») ® \A nv ) (2.6) 
k=lj=l 
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TABLE I: Reorderings of a matrix X representing an operator which acts on a composed Hilbert space. The arrows denote 
the indices exchanged. 



Transformation 


definition 


symbol 


preserves 
Hermicity 


preserves 
spectrum 


transposition 


Y T 

n.v 


— - X nv 
mu 


II 


yes 


yes 


flip 


X F 

nu 


— Xum 
vn 


^ — ^ 
<-> 


yes 


yes 


partial 


V^A 
nv 


— X nfi 
mv 


I- 


yes 


no 


transpositions 


Y Tb 

11.1/ 


Xmv 
nfi 


• i 


yes 


no 


reshuffling 




= Xmn 

flV 


/ 


no 


no 


reshuffling ' 


yR' 

nv 


— — X i/fj, 

n m 


\ 


no 


no 


partial 


X F 1 
nv 


= X tim 
nv 


<— > 


no 


no 


flips 


X F 2 
nv 


Xm.fi 

UTl 


<—> 


no 


no 



where Ckj = Tr((Ak <S> Aj)*X) may be neither Hermitian nor normal (which means that C and need not commute). 
The matrix X may be considered as a vector in the composite Hilbert-Schmidt space, 7Yhs ® T^hs, so applying its 
Schmidt decomposition |8l| we arrive at 

N 2 

|X)=£^K)®K'>, (2.7) 

where \f\k are the singular values of C, i.e. the square roots of the non- negative eigenvalues of CC^ . The sum of 
their squares is determined by the norm of the operator, Ylk=i ^ k = Tr(A"A^) = ||-X"||^ g . 

Since the Schmidt coefficients do not depend on the initial basis let us analyze the special case, in which the basis 
in Tins is generated by the identity matrix, U = 1 of size TV 2 . Then each of the TV 2 basis matrices of size N consist of 
only one non-zero element which equals unity, Ak — A mfl — \m){fi\, where k = N(m — 1) + /i. Their tensor products 
form an orthonormal basis in Tins ® 7~Lhs & n d allow to represent an arbitrary matrix X in the form l|2.6|l . In this 
case the matrix of the coefficients C has a particularly simple form, Cmp, — Tr[M mM ® A nv )X\ = Xmn. 

This particular reordering of a matrix deserves a name so we shall write X R = C(X) defining the following procedure 
of reshuffling of matrices. Using this notion our findings may be summarized in the following lemma: 

Schmidt coefficients of an operator X acting on a bi-partite Hilbert space are equal to the squared singular values 
of the reshuffled matrix, X R . 

More precisely, the Schmidt decomposition (|2.7|l of any operator X of size MN may be supplemented by a set of 
three equations 

{Afe}f=! = {SV{X R )} 2 : eigenvalues of (X R )^X R 

\A') : reshaped eigenvectors of [X R yX R , (2.8) 

\A") : reshaped eigenvectors of X R (X R )^ 

where we have assumed that TV < M. The initial basis is transformed by a local unitary transformation W a (£> Wb, 
where W a and Wb are matrices of eigenvectors of matrices (X R yX R and X R (X R y, respectively. Iff rank r of 
X R (X R y equals one, the operator can be factorized into a product form, X = X\ eg) Xi, where X\ = Tr2^T and 
X 2 = TnX. 

In general, one may reshuffle square matrices, if its size K is not prime. The symbol X R has a unique meaning if a 
concrete decomposition of the size K = MN is specified. If M ^ N the matrix X R is a TV 2 x M 2 rectangular matrix. 
Since (X R ) R = X we see that one may also reshuffle rectangular matrices, provid ed bo th dimensions are squares of 
natural numbers. Similar reorderings of matrices were considered by Hill et al. |80lll09| while investigating CP maps 
and later in pi llH 13^. I87L l88| to analyze separability of mixed quantum states. 

To get a better feeling of the reshuffling transformation observe that reshaping each row of an initially square 
matrix X of size MN according to Eq. (|2.1|) into a rectangular M x TV submatrix, and placing it according to the 
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lexicographical order block after block, one produces the reshuffled matrix X R . Let us illustrate this procedure for 
the simplest case N = M = 2, in which any row of the matrix X is reshaped into a 2 x 2 matrix 





X u 


X12 


X21 


X22 


x kj = 


X 13 


A14 


X23 


X 2 4 


X31 


X32 


-X4I 


-^42 




. X 33 


^34 


X43 


X44 _ 



(2.9) 



The operation of reshuffling could be defined in an alternative way, say the reshaping of the matrix A from l|2.1|) 
could be performed column after column into a vector a'. In the four indices notation introduced above (Roman 
indices running from 1 to N correspond to the first subsystem, Greek indices to the second one), both operations of 
reshuffling take the form 



A ,, 



Xmn 

^LV 



and 



X R n = Xi/f, 
nv nn 



(2.10) 



However, both reshuffled matrices are equivalent up to a certain permutation of rows and columns and transposition, 
so the singular values of X R and X R are equal. It is easy to see that (X ) R = X. In the symmetric case with M = N \ 
N 3 elements of X do not change their position during the operation of reshuffling (these are typeset boldface in 
Q2.9|) l: while the other N 4 — N 3 elements are exchanged. The space of complex matrices with the reshuffling symmetry, 
X = X R , is thus 2N 4 - 2{N 4 - N 3 ) = 2N 3 dimensional. 

For comparison we provide analogous formulae showing the action of partial transposition: with respect to the first 
subsystem, T4 = T ® 1 and with respect to the second, Tb = 1 <£> T, 



X m A n — Xnf, 

nv ml 



and 



v 1 b v 

nv n f J ' 



(2.11) 



Note that all these operations consist of exchanging a given pair of indices. However, while partial transposition (|2.11|l 
preserves Hermicity, the reshuffling l|2.10|l does not. For convenience we shall define a related transformation of flip 
among both subsystems, X^u = Xum, the action of which consists in relabeling of certain rows (and columns) of the 

nv vn 

matrix, so its spectrum remains preserved. Note that for a tensor product X = Y £g> Z one has X F — Z ®Y . In full 
analogy to partial transposition we use also two operations of partial Hip (see table QJ. All the above transformations 
are involutions, since performed twice they are equal to identity. It is not difficult to find relations between them, e.g. 
X p i = [{X R ') T *] R ' = [{X R ) T *] R . Since X R ' = [{X R ) F ] T = [(X R ) T } F , while X T * = (X T *) T and X F > = (X F *) F , 
thus the spectra and singular values of the reshuffled (partially transposed, partially flipped) matrices do not depend on 
the way, each operation has been performed, i.e. eig(X R ) = e\g(X R ) and SV(A' fl ) = SV(X R ), (eig(X Fl ) = eig(X F2 ) 
and SVpC^ 1 ) = SV(A" F2 )). 



III. COMPLETELY POSITIVE MAPS 



Let p € be a density operator acting on a TV-dimensional Hilbert space TLn- What conditions need to be 

fulfilled by a map $ : — > Ai^ N \ so it could represent a physical operation? If the map is linear than the image 
of a mixed state p does not depend on the way, how the state p was constructed out of projectors. This is a very 
important feature, since it allows for a probabilistic interpretation of any mixed state. Furthermore, we assume that 
the map is trace preserving, Trp = Tr^(p), which corresponds to the conservation of probability. 

Any quantum operation has to be positive - it should map any positive density operator into a positive operator. 
Strictly speaking the term 'positive' refers to positive semideRnite Hermitian operators, which do not have negative 
eigenvalues. However, this condition occurs not to be sufficient to produce physically realizable transformations. Any 
quantum state p may be extended by an ancilla cr into a tensor product acting in a KN dimensional Hilbert space. 
Hence the evolution of a linear map $ may be considered as the evolution of a part p of a larger system p eg) a. 
Therefore we should require complete positivity |9|, |90( , which means that for an arbitrary K dimensional extension 

Hn^'Hn®'Hk the map $®1 K is positive. (3.1) 

If the above requirement holds for any fixed K, such a property is called K-positivity [2i|, ||3(], |9(3 . Constructing the 
dynamical matrix we will show in the next section that any A^-positive map $ acting on Ai ( N ) is completely positive 
|18j . The importance of complete positivity in quantum mechanics was emphasized in the seventies by Kraus [59J, 
Lindblad |68( and Accardi and the restriction to CP maps is not trivial: There exist quantum maps which are 
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positive but not completely positive 90] . A simple and important example consists of transposition T of an operator 
in a given basis which for Hermitian operators is equivalent to complex conjugation: if A = then A T = A. 
Although the map T(p) = p T preserves the spectrum and therefore is positive, its extension Ta = T<g> 1, called partial 



transposition, is not, as discussed in section fVIII 

A linear CP-map which preserves the trace is called a quantum operation or quantum channel. Any linear, com- 
pletely positive map <!> may be represented [51^ . l60| by a collection of k Kraus operators Ai in the so-called Kraus 
form 



p^ p 



$(p)=J2A iP At 



(3.2) 



It is also called an operator sum representation of the map <E> or the Stinespring form, since its existence follows from 
the Stinespring dilation theorem |2y, |90| . If the set of Kraus operators satisfies the completeness relation 



k 

E 

z=l 



AlAi = 1, 



(3.3) 



the map is trace preserving: for any initial state p one has Tr(J2i=i A \ A iP) = Tr [ $ (/°)] = Tr P = 1 - Tmis . thc Kr aus 
operators may be considered as measurement operators while any trace preserving CP map may be interpreted as a 
generalized measurement [8l| . 

Let us denote the elements of the Kraus operators A; represented in an orthonormal basis by A 
matrix M of size N composed of non-negative entries [£ 



(i) 



. Define a Kraus 



i=l 



I An 



(3.4) 



where the Kraus operators are now denoted by A^ . The completeness relation (|3.3|l enforces that the sum of elements 
in each column of M equals to unity, so it is a stochastic matrix. Thus quantum operations are often called stochastic 
maps. 



system + environment 



d=N 2 k 2 -l 




b) 



environment 



partial tracing 
p , =Tr 2 (U(pg>c^)U , ) 



principal system 




principal system 



FIG. 1: Quantum operations represented by a) unitary operator U of size NK in an enlarged system including the environment, 
b) black box picture, in which U couples the principal ^-dimensional system p with the environment o of dimensionality K. 



Any map written in the Kraus form (|3.2|l is linear and preserves the trace if the condition (|3.3(l is fulfilled. Let p 
denote the density operator of the investigated quantum system. The time evolution of an isolated quantum system 
is unitary, while all sorts of a non unitary dynamics reflect its interaction with an environment, described by a density 
operator a . However, if we extend the system and study the dynamics of the total syst em, p ® a, composed of the 
system under investigation and the environment, its time evolution remains unitary [9tl27||. The non unitary operation 
of p emerges as an effect of the partial tracing with respect to the environment, 



p^p' = §{p) =Tr env [C7(p®a)C7t 



(3.5) 
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The above form is called as environmental representation of the map <I>. The entire process may be considered as 
a composition of three basic steps: adding an ancilla, unitary transformation, and partial tracing. The forms i|3.2[) 
and 1)3. 5J1 are equivalent in the sense that any quantum operation may be written either way. Both representations 
are not unique, but as discussed in section llVl the notion of a dynamical matrix allows one to define a distinguished, 
canonical Kraus form. 

Let the principal TV-dimensional system p be subjected to an operation defined by the kN dimensional matrix U 
and an initially pure state of the environment a = 

k 

p> = = Tr cnv [f/(p® \v){u\)rf] = J2(»\ U \v)p(»\ U ^)- ( 3 - 6 ) 

Since the fixed environmental state \v) belongs to the fc-dimensional Hilbert space 7i e nv, the expression {p\U\v) 
represents a square matrix of size TV, which we shall call A^. In the double index notation, introduced above, any 
element of the unitary matrix is denoted by J7m^, while the matrix An. — {p\U\u) consists of elements 

nu 

A%l = (m,u.\U\n,p) = U 7 n j f. (3.7) 

It is easy to see that 1)3.6)1 takes the operator sum form, Ylu=i A^pAj^. Moreover, due to unitarity of U the operators 
A^ satisfy the completeness relation 

k k 

J2 A l A » = E^VXmIE^) = (v\rfU\u) = In (3.8) 
i i 

and may be considered as Kraus operators. Every element of a given matrix A^ belongs to one of the TV 2 blocks of 
size k of the unitary matrix U. However, if we reorder it constructing U^ u — Ufim, then each matrix A„ is just a 

nv "tJ 

truncation of the flipped matrix U F ', since Amn = UF m represents its minor. 

un 

If the initial state of the environment in the representation 1)3. 5JI is chosen to have full rank, a = X^=i 9^l z/ )( l/ |) the 
operator sum representation consists of kN terms 

k kN 

v=l 1 = 1 

where Ai = ^/q^(p,\U\i'} and I = p + v{k — 1). 

In this way we have shown that for any operation written in the environmental form 1)3. 5J) we may find a corre- 
sponding Kraus form. Conversely, for any quantum operation in the Kraus form l|3.2|) consisting of k operators we 
may construct an environmental form |3, l67| . For instance take the first N columns of a unitary matrix U of size kN 
are constructed of k Kraus operators as defined in (|3.7() . Due to the completeness relation (|3.3|1 they are normalized 
and orthogonal. The matrix U has to be completed by N(k — 1) complex orthogonal vectors; they do not influence 
the quantum operation and may be selected arbitrarily. 

If the initial state of the environment is pure its dimensionality k of 7i cnv needs not to exceed N 2 , the maximal 
number of Kraus operators required. If the environment is initially in a mixed state, its weights coefficients qj are 
needed to specify the operation. Counting the number of parameters one could thus speculate that the action of any 
quantum operation may be simulated by a coupling with a mixed state of the environment of size N. However, this 
is n ot th e case: already for N = 2 there exist operations which have to be simulated with 3-dimensional environment 
|99l ll 1 1| , and the general question of the minimal size of Tt cnv remains open. 

It is also illuminating to discuss a special case of the problem, in which the initial state of the TV-dimensional 
environment is maximally mixed, a = p* = 1n/N. The unitary matrix U of size TV 2 , defining the map, may be 
treated as a vector in the composed Hilbert-Schmidt space Tins ® Hhs and represented in its Schmidt form (|2.7|l . 

U = X)i=i V^i\ A i) ® \ A 'i)i where Ai are eigenvalues of (U R yiI R . Since the operators obtained by reshaping 
eigenvectors of (U R yU R , form an orthonormal basis in Has, the procedure of partial tracing leads to the Kraus form 
consisting of TV 2 terms: 

p' ='i> u (p) = Tr cm [u(p®^l N )W 

N 2 N 2 N 2 

= Tr -v [E E [A lP A]) ® {jfAA'D] = - E \A iP A\. (3.10) 



i—1 j—1 i—1 
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Operations, for which there exist a unitary matrix U providing a representation in the above form, we shall call 
unistochastic channels. In a way, these maps are analogous to classical transformations given by unistochastic matrices, 
p 1 = Tp, where = |L/jj| 2 , since both dynamics are uniquely determined by a unitary matrix U - of sizes N and N 2 
in the classical and the quantum cases, respectively. 

In general one may consider maps analogous to (|3.10() with an arbitrary size of the environment. In particular we 
define generalized, K -unistochastic maps, determined by a unitary matrix U(N 1+K ), in which the environment of 
size N K is initially in the maximally mixed state, N~ k 1 N k. Such operations were analyzed in context of quantum 
information processing [54. l85| . and, under the name 'noisy maps', by studying reversible transformations from pure 
to mixed states |40) . By definition, 1-unistochastic maps are unistochastic. 

For any unistochastic map the standard Kraus form l|3.2[) is obtained by rescaling the operators, A; = \J\jjN At. 
Note that the matrix U is determined up to a local unitary matrix V of size N, in sense that U and U' — U(l ® V) 
generate the same unistochastic map, $£/ = $[/'■ 



IV. DYNAMICAL MATRIX 



A density matrix p of finite size N may be treated as a vector p reshaped according to (|2.1I) . The action of a linear 
superoperator <I> : p — > p 1 may thus be represented by a matrix L (" L" like linear) of size N 2 

p' = Lp or Pmu = L mt ip w , (4.1) 

nu 

where summation over repeated indices is understood. Nonhomogeneous linear maps p ' = L' p + a may also be 
described in this way. To obtain the homogeneous form (14.11) it suffices to substitute the matrix L' mu by Lmu = 

nu nv 

nu 

We require that the image p' is a density matrix, so it is Hermitian, positive, and normalized. These three conditions 
impose constraints on the matrix L: 



i) P ' = (p') t =► L tl m = L* mli so L* = L F , (4.2) 

ii) p 1 > ==>■ Lmup nu > for any state p, (4.3) 



N 



111 



trp' = 1 => ^2 L ™™ = S w ( 4 - 4 ) 



Note that l|4.2|l is not the condition of Hcrmicity and in general the matrix L representing the operation $ is not 
Hermitian. However, if we reshuffle it according to l|2.10|l and define the dynamical matrix 

D$ = L R so that Dmu, = Lmn, (4.5) 

nv M" 

than Dg, is Hermitian, D = D' , due to (|4.2|l . The linear superoperator L is then uniquely determined by the dynamical 
matrix, since L = D R . The notion of dynamical matrix was introduced already in 1961 by Sudarshan, Mathews and 
Rau [94| . Later such matrices were used by Choi [n| and are sometimes called Choi matrices j3g . 

What conditions must be satisfied by a Hermitian matrix D of size N 2 to be a dynamical matrix? The trace 
condition (|4.4H . rewritten below for the matrix D, determines its trace 

N N N N 

Tr A D = V Dmu = 5^ TrD = V V D m u = V = N. (4.6) 

— * mv — * — * rrifi — * 

The positivity condition Ij4.3fl implies that for any states \x) and p the expectation value (x\D R p\x) is not negative. 
Assuming that the initial state is pure, p = \y)(y\ so that p^ = y^y*, we obtain 

(x\D R (\y)(y\)\x) = ((x\ ® (y\)D(\y) ® |x» > 0, (4.7) 

since the double sum DmuyuVt represents as well the matrix D R (\y)(y\) as well as (y\D\y). Thus the dynamical 

nu 

matrix D must be positive on product states \y) <g> \x). This property is called block-positivity. In fact Jamiolkowski 
proved, in 1972, that the converse is also true: if property (|4.7|l is satisfied then the map $ determined by matrix D 
is positive - see section IVIIII 
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Positivity of p$> is a sufficient, albeit not necessary requirement for i|4.7[l . If £>$ > then the map <E> is completely 
positive 0, |2jJ [3g, 0, US • To prove this let us represent by its spectral decomposition, 

k k 

D*=J2 d ^ (*< I 80 that ^ = E ^X^xS ■ (4-8) 

i=l i=l 

Here fc denotes the rank of D$ so k < N 2 , while ^ represent N x N matrices obtained by reshaping the eigenvectors 
of length TV 2 . Due to (|4.6I) the sum of all eigenvalues di is equal to N. If Z?$ is positive, all its eigenvalues are 
nonnegative, so we may rescale the eigenvectors defining the operators A — A 1 ) 

A%=^d lX %, i=l,2,...,fc, (4.9) 

so that 

k k 

D W = £ A % A nl = £(^ ® ( 4 -!0) 

i=l i=l 

If all Kraus operators are real A = A, or purely imaginary, A; = — A;, then the Hcrmitian dynamical matrix is real 
and hence symmetric, D<& = D^. Reshuffling the dynamical matrix D (of the size N 2 , with a finite N), we may write 



k k 

L^i.d,^^^ 1 ' 1 . (4-11) 

i=l i=l 

The latter form may be thus considered as a Schmidt decomposition (|2.7() of the superoperator L for a CP-map, since 
the non- negative eigenvalues of D are simultaneously Schmidt coefficients of L = D R . 

Alternatively one may find a matrix A of size ./V 2 x k such that D$ = AA^. This is always possible since D is 
positive. Then the operator A z ) is obtained by reshaping the i-th column of A into a square matrix. Using the eigen 
representation of Z?$ the quantum map $ defined in (|4.1|) becomes 

k 

p'm* = Drnnp nv = £ A^ nPnv A% (4.12) 
i=l 

and it may be written in the canonical Kraus form 

k k 

^ = X> x {i) p(x {t) V = £>/>4 (4 - 13) 

i=l i—1 

Let us compute the matrix E = y^-, A\Ai, 

k k 

E,u = £ A« A» = £ dkx%X { l = Drn,, (4.14) 

2=1 2=1 

and therefore E = Tr^Z?^. Due to the trace preservation constraint l|4.6|) E = 1, so the operators A satisfy the 
completeness relation (|3.3|l . Hence l|4.13|l is equivalent to the Kraus form (|3.2|l and represents a trace preserving CP 
map. We have thus shown that any positive dynamical matrix D, which satisfies the condition (|4.6|l . specifies uniquely 
a quantum operation, Hence any iV-positive map $ leads to a positive matrix D$, so $ is completely positive. 

On the other hand, the Kraus representation l|3.2|l is not unique - two sets of the Kraus operators Ai,i = 1, 
and Bj,j — 1, n represent the same operation if and only if the dynamical matrices given by l|4.10|l are equal. This 
is the case if there exists a unitary matrix V of size mN such that A = BV so that the dynamical matrices both sets 
generate, are equal, 

D = AA f = BV(BV) f = BB f (4.15) 

Here m = max{?,n} and the shorter list of the Kraus operators is formally extended by \l — n\ zero operators. 

In principle the number of Kraus operators in l|3.2|l may be arbitrarily large. However, for any operation one may 
find its Kraus form consisting of not more operators than the rank k of the dynamical matrix D. This Kraus rank 
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TABLE II: Quantum operations $ : MS ' —* MS . properties of superoperator L and dynamical matrix = L 



Matrices 


Superoperator L — D R 


Dynamical matrix D$> 


Hermicity 


No 


Yes 


Trace 


spectrum is symmetric 
=> trL G R 


trD& = N 


(right) 
Eigenvectors 


invariant states 
or transient corrections 


Kraus operators 


Eigenvalues 


\Zi\ < 1 
— In - decay rates 


weights of Kraus 
operators, di > 


Unitary evolution 
Dif, = {U®U*) R 


||L|| 2 =JV 


S(<?) = 


Coarse graining 


|L|| 2 = VN 


S(d') = In AT 


Complete depola- 
risation, Dj. = 1 


= 1 


S{d') = 2 In TV 



will never exceed the dimension of the dynamical matrix equal to N 2 . Moreover, the Kraus operators Ai may be 
chosen to be orthogonal Q, |3(| . To find such a canonical Kraus form given by (|4.13|l for an arbitrary operation <f> it 
is enough to find its linear matrix L, reshuffle it to obtain the dynamical matrix diagonalize it, and out of its 
eigenvalues di and reshaped eigenvectors \xi) construct by (|4.9(l the orthogonal Kraus operators A^. They satisfy 

(A t \A 3 ) = TrA^Aj - ^d~d~( Xl \ X] ) = diSij. (4.16) 

In a sense the Kraus form 13.2fl of a quantum operation <I> can be compared with an arbitrary decomposition of 
a density matrix, p = ^Pi|^t)(^»|) while its eigen decomposition corresponds to the canonical Kraus form of the 
map, for which ||A,-|| 2 = dj. In the generic case of a nondegenerate dynamical matrix D, the canonical Kraus form is 
specified uniquely, up to free phases, which may be put in front of each Kraus operator Ai. 

As discussed in section IVIIII such an analogy between the quantum maps (dynamics) and the quantum states 
(kinematics) may be pursued much further. Let us state at this point that the dynamical matrix is a linear function 
of the quantum operations $ and $ in the sense that 

D a$+fc * = aL>$ + bD^. (4.17) 

An arbitrary quantum operation $ is uniquely characterized by any of the two matrices L or _D$ = L R , but the 
meaning of their spectra is entirely different. The dynamical matrix £)$ is Hermitian, while L is not and in general 
its eigenvalues Zi are complex. Let us order them according to their moduli, \z\\ > \z-z\ > • • • I^a^I > 0. A continuous 
linear operation $ sends the convex compact set MS N ' into itself. Therefore, due to the fixed-point theorem, this 
transformation has a fixed point - an invariant state u\ such that Lu\ = <j\ . Th us Z\ = 1 and all eigenvalues fulfill 
\z%\ < 1, since otherwise the assumption that $ is positive would be violated |l02j . 

The trace preserving condition applied to the eigen equation Lai = Zi<Ji, implies that if Zi ^ 1 than tr(<7;) = 0. If 
r = \ z 2\ < 1 then the matrix L is primitive jT^ and all states converge to the invariant state o\. If L is diagonalizable 
(there is no degeneracy in the spectrum or there exists no nontrivial blocks in the Jordan decomposition of L, so that 
the number of right eigenvectors (Tj is equal to the size of the matrix N 2 ), then any initial state po may be expanded 
in the eigenbasis of L, 

JV 2 N 2 

Po = ^ CjQ-j while p t = L'po = ^CjZ-cTi. (4-18) 

i=l i=l 

Therefore po converges exponentially fast to the invariant state <j\ with the decay rate not smaller than — In r and 
the right eigenstates <7j for i > 2 play the role of the transient traceless corrections to po- The super-operator L sends 
Hermitian density matrices into Hermitian density matrices, p\ = p\ = Lpo = Lp , so 

if Lx = Z X then Lx* = z*x\ (4.19) 

and the spectrum of L (contained in the unit circle) is symmetric with respect to the real axis. Thus the trace of L 
is real, as follows also from the Hermicity of Dg, = L R . 
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On the other hand, the real eigenvalues dj of the dynamical matrix D$, satisfy the normalization condition (14. 6|) . 
dj ~ ^> wm ch we assume to be finite. If the map $ is completely positive all eigenvalues dj of D$ are non- 
negative, and the matrix p$ = may be interpreted as a density matrix acting in 7i N 2. The eigenvalues of p$, 
equal to d' = d/N, determine the weights of different operators Aj contributing to the canonical Kraus form of the 
map given by l|4.9|l . To characterize this probability vector quantitatively we use the Shannon entropy to define the 
entropy of an operation $, 

N 2 

S(9) = S(—d) = - Y,( d i/ N ) Hdi/N) = S N {p$), (4.20) 

i=l 

equal to the von Neumann entropy of p$. In order to characterize, to what extent the distribution of the elements of 
the vector d is uniform one may also use other quantities like linear entropy, participation ratio or generalized Renyi 
entropies. If S(d') — 0, the dynamical matrix D$ is of rank one, so the state pg> is pure. Then Eq. (|4.10|l reduces 
to Dmii = U m JJ*, while L — U ®U* and the map <I> represents a unitary rotation. The larger the entropy 5* of an 
operation, the more terms enter effectively into the canonical Kraus form, and the larger are effects of decoherence in 
the system. 

For any finite N the entropy is bounded by S max = 21niV, which is achieved for the rescaled identity matrix, D = 
t N 2/N. This matrix represents the completely depolarizing channel $*, such that any initial state p is transformed 
into the maximally mixed state, 

$*(p) = P* = 1n/N for any state p. (4-21) 

Under the action of this map complete decoherence takes place already after the first iteration. It is easy to see that 
the spectrum of the corresponding superoperator L* = jj(1n) r consists of Z\ = 1 and Zj = for i > 2. In general 

the norm of the superoperator, ||L||2 = vTrLlJ, may thus be considered as another quantity characterizing the 
decoherence induced by the map. The norm ||£||2 varies from unity for the completely depolarizing channel <&» (total 
decoherence) to N for any unitary operation (no decoherence). Both quantities are connected, since equal to 

||-D||| = trZ? 2 , is a function of the Renyi entropy of order a = 2 of the spectrum of D rescaled by 1/N. Making use 
of the monotonicity of the Renyi entropies ^1 we S e * ll-^lk > iVexp( — 5(<I>)/2). 

An alternative way to characterize the properties of a superoperator is to use its trace norm, 

N 2 

7 = ||L|| X =Tr|L| =TtVlU=J2^ ( 4 - 22 ) 

i=i 

where are the singular values of L = D R . As discussed in this quantity is useful to determine separability 

of the state associated with the dynamical matrix, p$ = = D/N. 

If the set of Kraus operators determines a quantum operation $ then for any two unitary matrices V and 

W of size N the set of operators A[ = VAiW satisfies the relation (|3.3|) and defines the operation 

k k 

P^p' = $vw(p) =£4/^ = v(J2MWpW^)Al)vl (4.23) 

i=l i=l 

The operations and <&vw are in general different, but unitarily similar, in a sense that the spectra of the dynamical 
matrices are equal. Thus their generalized entropies are equal, so S(Q) = S($vw) an d \ \L\U = \ \Lvw\\2- The latter 
equality follows directly from the following law of the transformation of superoperators, 

L vw = (V® V*)L(W® W*), (4.24) 

which is a consequence of l|4.11|l . 



V. UNITAL & BISTOCHASTIC MAPS 



Consider a completely positive quantum operation $ defined by Kraus operators l|3.2|) . If the set of k operators Ai 
satisfies the relation, 



k 



z ^A l A\ = \ = tx B B (5.1) 



11 



the operation is united (or exhaustive), i.e. the maximally mixed state p* = I/N remains invariant, $(p*) = 
w S,'— 1 Ai-Af = p*. Since (p*) mn = iiJ mn then the elements of $(p*) are -hLmnSnu — TfDmn, which explains 
the right hand side equality in i|5.1f) . related to the properties of the dynamical matrix .D. 

Observe the similarity between the condition (|3.3|l for the preservation of trace and the unitality constraint (|5.1|l . 
Unitality imposes that the sum of all elements in each row of the corresponding matrix M defined by (|3.4jl is equal 
to one. Hence for any map $ which is simultaneously trace preserving and unital, the Kraus matrix M is doubly- 
stochastic (bistochastic) . Therefore a trace preserving, unital completely positive map is called a bistochastic map 
[^.l64j. and may be considered as a noncommutative analogue of the action of a bistochastic matrix Bona probability 
vector - see table Hill In the former case, the maximally mixed state p* is ^-invariant, while the uniform probability 
vector is an invariant vector of B. 

If all Kraus operators are Hermitian, Ai — A\, the channel is bistochastic but this is not a necessary condition. 
For example, unitary evolution may be considered as the simplest case of the bistochastic map with k = 1. A more 
general class of bistochastic channels is given by a convex combination of unitary operations, also called random 
external Gelds (REF) |(J, 

fe fe 
ft = ®{p) = Y J PiV i pV} , with Pi >0 and 5^ = 1, (5.2) 

where each operator V$ is unitary. The Kraus form 13. 2|) can be reproduced setting Ai = y/plVi- 

In general, a map for which all Kraus operators are normal, [Ai, A\\ = 0, is bistochastic. Any convex combination 
of two bistochastic maps is bistochastic, and similarly, any convex combination of two random external fields belongs 
to this class. Thus the spaces of bistochastic maps and the spaces of random external fields are convex. For N = 2 
both sets coincide and any bistochastic map may be represented as a combination of unitary operations - see Fig. [^i. 
Such one-qubit bistochastic maps are also called Pauli channels. For N > 3 the set of bistochastic maps Bn is larger 
than the set of REFs 0. 

For any quantum channel $(p) one defines its dual channel $(p), such that the Hilbert-Schmidt scalar product 
satisfies (<&(cr)|p) = (cr|$(p)) for any states a and p. If a CP map is given by the Kraus form <E>(p) = J2i A%pA\, 
the dual channel reads &(p) = Y^i A\pAi. Making use of (|4.10|) we obtain a link between the dynamical matrices 
representing dual channels, 

L 4 = (Llf = (iff and = (D%) F = (D F f = D$ . (5.3) 

Since neither the transposition nor the flip modify the spectrum of a matrix, the spectra of the dynamical matrices 
for dual channels are the same, as well as their entropies, S(Q) — S(&). 

Comparing the conditions l|3.3|l and (|5.1I) we see that if channel $ is trace preserving, its dual <t is unital, and 
conversely, if channel $ is unital then <& is trace preserving. Thus the channel dual to a bistochastic one is bistochastic. 

For any quantum operation $ given by the Kraus form l|3.2|) one defines its effect by E = X)j=i AiA\. Due to Eq. 
(|5.1|) the effect is obtained by partial trace of the dynamical matrix, E = TtbD, so it is Hermitian and positive. Since 
E(Q) — jf&(p*), it is clear that the effect of any operation does not depend of the particular choice of the Kraus 
operators used to represent it. For bistochastic maps 3>(p*) = p* and E = 1. 

Bistochastic channels are the only ones which do not decrease the von Neumann entropy of any state they act on. 
To see this consider the image of the maximally mixed state, with maximal entropy In N. If the map is not unital, 
(the channel is not bistochastic), then p' = $(/?*) differs from p*, so its entropy decreases. As an important example 
of bistochastic channels consider the coarse graining operation, which sets all off-diagonal elements of a density matrix 
to zero, ^gg(p) = diag(p). It is described by the diagonal dynamical matrix, .D~~ = SmnS^Sm^, with N elements 

nv 

equal to unity. 

Let us analyze in some detail the set BIA n of unistochastic operations, for which there exists the representation 
(|3.1U|I . The initial state of the environment is maximally mixed, a — p* = 1/N, so the quantum map tyy is determined 
by a unitary matrix U of size TV 2 . The reshaped Kraus operators Ai are proportional to the eigenvectors 14.9(1 of the 
dynamical matrix D^, u . On the other hand, they enter also the Schmidt decomposition l|2.7l) of U as shown in l|3.10|l . 
and are proportional to the eigenvectors of (U R yU R . Therefore 

D^ = ^(UyU R so that L = ±[(U R yU R ] R . (5.4) 

We have thus arrived at an important result: for any unistochastic map the spectrum of the dynamical matrix D 
is given by the Schmidt coefficients, dj = Xi/N, of the unitary matrix U treated as an element of the composite HS 
space. Hence the entropy of the operation S(^u) is equal to the entanglement entropy of the unitary matrix, S(U). 
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TABLE III: Quantum maps acting on density matrices and given by positive definite dynamical matrix D > versus classical 
Markov dynamics on probability vectors defined by transition matrix T with non-negative elements 



Quantum 


completely 
positive maps: 


Classical 


Markov chains 
given by: 


S? 


Trace preserving, TrA-D = 1 


S? 1 


Stochastic matrices T 


s§ 


Unital, Tr s L> = 1 




T T is stochastic 


cQ 


Unital & trace preserving maps 


^3 


Bistochastic matrices B 


s? 


Maps with A4 = A\ 
, , ■ =? T)= P T . 




Symmetric bistochastic 
matrices, B — B T 


sf 


Unistochastic operations, 
D = U R (U R ) j 




Unistochastic matrices, 

Bij — \Uij\ 2 


„q Orthostochastic operations, 
*6 D = qR {0 R )T 


s 6 Ci 


Orthostochastic matrices, 
B^ = \Oij 2 


s? 


Unitary transformations 


s? 1 


Permutations 



Moreover, the linear entropy of entanglement of U studied in |ll2j| is a function of the norm of the superoperator, 
E(U) = 1 — J2i = 1 — ll-^lli/^ 2 - ^ vanishes for any local operation, U = XJ\ ® U2, for which the superoperator 
is unitary, L — U\ ®Ul so \\L\\ 2 — N 2 . The resulting unitary operation is an isometry, and can be compared with 
a permutation S7 acting on the simplex Ajv-i of classical probability vectors. If the matrix U is orthogonal the 
corresponding dynamical matrix is symmetric, D^ = {O r ) t O r = D^ o . The corresponding operation will be called 
an orthostochastic map, as listed in Table ITTT1 The spaces listed there satisfy the following relations 5ifl52 = S3 and 
S3 D S5 D Sg and S5 D S7 in both, classical and quantum set ups. However, the analogy is not exact: the inclusions 
S^ 1 C S3 1 and S? 1 C Sq 1 and S$ C S® do not have their counterparts. Note that unitary (orthogonal) matrices 
defining quantum maps M.( N > — > are of size N 2 while these determining Markov chains Aat_i — > A^-i are of 

size N. Note that already for N = 2 not all bistochastic maps are unistochastic [73j . 



VI. ONE QUBIT MAPS 



In the simplest case N = 2 the quantum operations are called binary channels. In general, the space CPn is 
N 4 — N 2 dimensional. Hence the space CP2 of binary channels has 12 dimensions and there exist several ways to 
parametrise it |29t l53j . The dynamical matrix provides a straightforward, but not very transparent method to achieve 
this goal. Any Hermitian matrix D of size 4, which satisfies tr aD = 1 may be written as 



D = 



1 

2 


+ a 




X 




y 


z 




X 


1 

2 


- b 




w 


—u 




y 




w 


1 

2 


— a 


—x 




z 




-u 




—x 


h + K 



(6.1) 



where a and b are real and x,y, z,u,w complex parameters. If they are adjusted in such a way to assure positivity 
of D, then this matrix represents a trace preserving CP map. Note that D/N represents a certain density matrix in 
M^ 4 \ If additionally a = b and y — u then the condition tr^D = 1 is fulfilled, so the map is bistochastic. 

An alternative approach to the problem is obtained using the Stokes parametrisation, which involves Pauli matrices 
Ofc. Any state p = (1 + f ■ <?)/2, is characterized by the Bloch vector f. If p' = Q(p) is represented by f', then any 
linear dynamics inside the Bloch ball may be described by an affine transformation of the Bloch vector 



t = tf + K, 



(6.2) 



where t denotes a real matrix of size 3, while k = (n x , n y , k z ) is the translation vector. It determines the image of 
the ball center, ^(/c*) = (1 + K ■ er)/2. For unital maps &(p*) = p* so k = 0. Any positive map sends the Bloch 
ball into an ellipsoid which may degenerate to an ellipse, an interval or a point. Out of 12 real parameters in l|6.2|) 
3 eigenvalues of tv determine the size of the ellipsoid, 6 related to the eigenvectors specify the orientation of its axis 
and the 3 parameters of the vector k characterize the position of its center. 

The real matrix t may be brought to the form, t — 0\r\Oi, where Oi € 0(3) and the diagonal matrix 77 contains 
non-negative singular values of t. However, it is convenient to impose an extra condition that both orthogonal matrices 
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represent proper rotations (Oj £ 5*0(3)), at the expense of allowing some components of the diagonal vector fj to 
be negative. This decomposition is not unique: the moduli \r]i\ are equal to singular values of t and the sign of the 
product r) x r) y ri z = det(t) is fixed. However, the signs of any two components of fj may be changed by conjugating the 
map with a Pauli matrix. 

In general, for any one-qubit CP map $ one may find unitary matrices U and W of size 2 such that the unitarily 
similar operation §uw defined by (|4.23l) is represented by the diagonal matrix rj. For simplicity we may thus restrict 
ourselves to the maps for which the matrix t is diagonal and consists of the damping (distortion) vector, fj = (r] x , T) y ,r) z ). 
In this case the ellipsoid has the form 



rx-n x y r y-K y y / z-k z \- 



= 1. 



(6.3) 



The affine transformation l|6.2|l determines the superoperator L of the map. Reshuffling the superoperator matrix 
according to l|4.5|) we obtain the dynamical matrix D = L R which corresponds to the map 



D 



- + Vz 




o 



y 



1 - Vz 
t]x - 



K, 



Vy 

IK 



1 



Tjx - Vy 
— T) z -K 



Vx + Vy 



1 + l] z - K z 



(6.4) 



f]x ~\~ f]y K x IKy 

clearly a special case of Ij6.1|l . For unital maps k — and the matrix D splits into two blocks and its eigenvalues are 

e£o,3 = g I 1 + Vz ± (Vx + Vy)] and d i,2 = ^[1 -Vx ± {Vx ~ Vv)]- ( 6 - 5 ) 
Hence if the Fujiwara-Algoet conditions [2|j 

(l±Vz) 2 >(v x ±Vy) 2 (6.6) 

are fulfilled, the dynamical matrix D is positive definite and the corresponding positive map $^ is CP. This condition 
shows that not all ellipsoids located inside the Bloc h ball may be o btained by acting on the ball with a CP map - see 
recent papers on one qubit maps pi IH1 l89l logt llM ITol fl06| . 

Note that dynamical matrices of unital maps do commute, D$J — 0. Hence they share the same set of 

eigenvectors. When reshaped they form the orthogonal set of Kraus operators consisting of the identity cr = la> an d 
the three Pauli matrices a. Thus the canonical form of an arbitrary one-qubit bistochastic map reads 



i=0 



(6.7) 



which explains the name Pauli channels. 

For concreteness let us distinguish some one-qubit channels. We are going to specify their translation and distor- 
tion vectors, k and fj, since with use of the transformation l|4.23|l we may bring the matrix t to its diagonal form. 
Alternatively, the channels may be defined using the canonical Kraus form l|4.13|l in which Kraus operators are given 
by the eigenvalues and eigenvectors of the dynamical matrix, Aj = ^/dixi ■ From l|6.7|l it follows that for any N = 2 
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TABLE IV: Exemplary one-qubit channels: distortion vector ff, translation vector k equal to zero for unital channels, rescaled 
Kraus spectrum d' , and the Kraus rank k. 



Channels 


ff 


K 


unital 


«? 


k 


rotation 


(1,1,1) 


(0,0,0) 


yes 


(1,0,0,0) 


1 


phase flip 


(l-p,l-p,l) 


(0,0,0) 


yes 


(l-p/2,p/2,0,0) 


2 


decaying 


(V1-p,V1-p,i-p) 


(0,0, p) 


no 


(l-p/2,p/2,0,0) 


2 


depolarizing 


[1-x] (1,1,1) 


(0,0,0) 


yes 


j(4 — 3x, x, x, x) 


4 


linear 


(0,0, q) 


(0,0,0) 


yes 


i(l + q,l-q, 
l-q,l+q) 


4 


planar 


(0,r,q) 


(0,0,0) 


yes 


j{l + q + r,l-q-r, 
1 — q + r, 1 + q — r) 


4 



bistochastic map the Kraus operators are Ai = \fd^ai. For the decaying channel, (also called amplitude-damping 
channel), which is not bistochastic, the eigenvectors of the dynamical matrix give 



A 1 = 







1 



V^~P 



and Ao — 



o Vp 





(6.8) 



Basic properties of some selected maps are collected in Table llVl and illustrated in Fig. [21 



The geometry of the set of bistochastic maps is simplest to understand for N — 2. As discussed in section Ivl this 
set coincides j£4( with the set of N — 2 random external fields l|5. 2JI , which may be defined as the convex hull of four 
unitary operations, (Jq = l,a x ,a y and o z . Thus, in agreement with l|6.7|l . any one-qubit bistochastic map <£> € B2 may 
be written as a convex combination of four matrices aj, j — 0, 3. Since Uj — — i exp(i7rtTj/2), and the overall phase 
— 7r is not relevant, each Pauli matrix represents the rotation of the Bloch ball around the corresponding axis by angle n. 
The distortion vectors if of four extremal maps <jj read (1, 1, 1), (1, —1, —1), (—1, 1, —1), and (—1, —1, 1), respectively. 
Hence, in accordance with the constraints Ijtj.bfl . the set of bistochastic maps B2 written in the canonical form (with 
diagonal matrix t), forms a tetrahedron, see Fig. Its center is occupied by the completely depolarizing channel 
with ff = 0, which may be expressed as a uniform mixture of four extremal unitarities; 'J* = (1 + a x + cr y + a z )/A. 
Interestingly, the set of non unital channels with a fixed translation vector k ^ forms a convex set which resembles 
a tetrahedron with its corners rounded H3 . 




,{1,1, 1) 0-1,1) 

T C ) tr— -nT 



( 1-1,1 



(-1,1,-1) 



(-1,-1 



(i,i^n_ 



1) 



-1,-1) 



n 



(-1-1-1) 



FIG. 3: Subsets of one-qubit maps: a) set B2 of bistochastic maps (unital & CP), b) set T{B2) of unital & CcP maps, c) set 
of positive (decomposable) unital maps. 



Observe that in the generic case \r]i\ < 1, which suggests that a typical one-qubit CP map is a contraction. This is 
true also for higher dimensions, N > 2. The monotone distances (e.g. the trace distance and the Bures distance |74j1 
do not increase under the action of CP maps. Apart from unitary (antiunitary) operations, for which these distances 
are preserved, the transformation inverse to a quantum operation is not an operation any more: some mixed states 
are sent outside the set of positive operators. 
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VII. POSITIVE & DECOMPOSABLE MAPS 

Quantum transformations which describe physical processes arc represented by completely positive maps. Why 
should we care about maps which are not completely positive? On one hand it is instructive to realize that seemingly 
innocent transformations are not CP, and thus do not correspond to any physical process. On the other hand maps 
which are positive, not completely positive provide a crucial tool in the investigation of entangled mixed states 

Consider the transposition of a density matrix in a certain basis T : p — > p T . The corresponding superoperator 
Lt entering (14.111 has the form (Lt)tom = 3 ml/ 5 n n, and equals the dynamical matrix, Lt = L% — Dt- This is a 
permutation matrix which contains N diagonal entries equal to unity and N(N — l)/2 blocks of size two, Thus its 
spectrum, spec(Z>r), consists of N(N + l)/2 eigenvalues equal to unity and N(N — l)/2 eigenvalues equal to — 1, 
which is consistent with the constraint trD — N. The matrix Df is not positive, so the transposition T is not 
completely positive. Another way to reach this conclusion is to act with the extended map of partial transposition on 
the maximally entangled state Q8.2[l and to check that [T <S> (ip\) has negative eigenvalues. 

The transposition of an iV-dimensional Hcrmitian matrix D, which changes the signs of the imaginary part of the 
elements Dij, may be viewed as a reflection with respect to the N(N + l)/2 — 1 dimensional hyperplane, so the 
remaining N(N — l)/2 components do change their signs. As shown in Fig. |I]this geometrical interpretation is simple 
to visualize for N = 2: the transposition corresponds to a reflection of the Bloch ball with respect to the (x, z) plane - 
the y coordinate changes its sign. Note that rotation of the Bloch ball around the z-axis by the angle 7T, realized by a 
unitary operation, p — > a z pa z , also exchanges the 'western' and the 'eastern' hemispheres, but is completely positive. 

z A 




FIG. 4: Non-contracting transformations of the Bloch ball: a) transposition (reflection with respect to the x — z plane) - not 
completely positive; b) rotation by 7r around z axis - completely positive. 

As discussed in section Hvl a given map $ is not CP, if the corresponding dynamical matrix D contains a negative 
eigenvalue. Let m > 1 denote the number of the negative eigenvalues (in short, the neg rank of D). Then the spectral 
decomposition of D takes the form 

N 2 -m N 2 

D= *IXi)(Xi| - E d i\Xi)(Xi\ ■ (7-1) 

i=l i=N 2 -m+l 

In analogy to the Kraus form (|4.13|l we may write the canonical form of a not completely positive map 

N 2 -rn N 2 

p'= £ ^X (i) p(x W ) f - E diX®P<X {i) ?, (7-2) 

i=l i=N 2 -m+l 

where the Kraus operators Ai = \fdix^ form an orthonormal basis. The above form suggests that a positive map may 
be represented as a difference of two completely positive maps 95] . Even though this statement is correct (for finite 
iV"), it does not solve the entire problem: taking any two CP maps and constructing a quasi- mixture (with negative 
weights allowed) , $ = ( 1 + a) #f p — P j we do not know in advance how large the contribution a of the negative par t 
might be, to keep the map $ positive... Although some criteria for positivity are known for several years |49l l7Ct l91 | . 
they do not allow one to perform a practical test, whether a given map is positive. A recently proposed technique of 
extending the system (and the map) a certain number of times is shown to give a constructive test for positivity for a 
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large class of maps 24] . In fact the characterization of the set Vn of positive maps: — > M.( N ^ is by far not simpl e 
and remains a subject of considerable mathematical interest [3 E3, E^, EE Hj, H3, E3, [Til IzS HH IlH HE EM ■ 

By definition, Vn contains the set CVn of all CP maps as its proper subset. 

To learn more about the set of positive maps we will need some other features of the operation of transposition T. 
For any operation <I> the modifications of the dynamical matrix induced by a composition with T may be described 
by the partial transpose transformation 

L T * = L F \ D T z = Dl A , and L $T = , D^ T = D% B . (7.3) 

To demonstrate this it is enough to use the explicit form of Lt and the following law of composition of dynamical 
matrices. Since the composition of two maps \&<i> results in the product of linear matrices, 

= L^Lq,, hence = [D§ D R ] R . (7.4) 

Even though the composition of two maps is usually written as W • $, to simplify the notation the symbol • will often 
be dropped. Note that both operations commute, = if [D R , D R ] — 0. 

Positivity of follows also from the fact that the composition of two CP maps is completely positive. Alterna- 
tively one may prove the following reshuffling lemma 

Consider two Hermitian matrices A and B of the same size KN. 

If A>0 and B > then (A R B R ) R > 0. (7.5) 

It was formulated in a different set up and proved in [3(| . 

The sandwiching of $ between two actions of transpositions does not influence the spectrum of the dynamical 
matrix, Lt$t = L$ = L*^ and Dt$t = D$ = D%. Thus if $ is a CP map, so is TQT, (if Z)$ is positive so is D]^) - 
see Fig. |3| 




FIG. 5: a) Sketch of the set of CP maps and of its image with respect to transposition: the set of CcP maps, b) the isomorphic 
sets of quantum states (dynamical matrices) and its image under the action of partial transposition. 



The non completely positive map T of transposition allows one to introduce the following definition |19|, |2lJ, |9l| : 
A map $ is called completely co-positive (CcP), if the map T<I> is CP. 

Properties l|7.3|l of the dynamical matrix imply that the map <I>T could be used instead to define the same set of 
CcP maps. Thus any CcP map may be written in a Kraus-like form 

p' = Y,A iP T A\. (7.6) 

Moreover, as shown in Fig. [S] the set CcP may be understood as the image of CV with respect to the transposition. 
Since we have already identified the transposition with a kind of reflection, it is rather intuitive to observe that the 
set CcP is a twin copy of CV with the same shape and volume. This property is easiest to analyze for the set B2 of 
one qubit bistochastic maps |7^|. The dual set of CcP unital one qubit maps, T(Bq) 1 forms a tetrahedron spanned by 
four maps Tai for i = 0,1, 2, 3, which is the reflection of the set of the bistochastic maps with respect to the center of 
the tetrahedron - the completely depolarizing channel <&* -see Fig. |3|d. Observe that the corners of B2 are formed by 
proper rotations, for which det(i) = ?7i?72?73 is equal to +1, while the extremal points of the set of CcP maps represent 
reflections for which det(t) = — 1. 
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A positive map $ is called decomposable, if it may be expressed as a convex combination of a CP map and a CcP 
map, $ = a$cp + (1 — a )$CcP with se [0,1]. A relation between CP maps acting on quaternion matrices and the 
decomposable maps defined on complex matrices was shown by Kossakowski An important characterization of 
the set V 2 of positive maps acting on (complex) states of one qubit follows from the St0rmer-Woronowicz theorem 

HJEil 

Every one-qubit positive map \& G V2 is decomposable. 

In other words, the entire set of one qubit positive maps can be represented by the convex hull of the set of CP and 
CcP maps, V2 = conv hull (CV2 U CcV^^j ■ This property, illustrated in Fig. holds also for the maps MS^ — > A^ 3 -* 

and — > M. 1 ^ |l07j . but is not true in higher dimensions, in particular in the sets Vn with N > 3. Consider a 
map defined on depending on three non-negative parameters, 



apu + bp 2 2 + cp 33 

Cpxi + ap22 + bp33 

bp u 






Cp 2 2 + aP33 



ft- 



(7.7) 



The map ^2, 0.2 € V3 was a first example of a indecomposable map found by Choi in 1975 |20| in connection with 
positive biquadratic forms. As denoted schematically in Fig. this map is extremal and belongs to the boundary of 
the convex set V3. The Choi map was generalized later in [22] and in [ljj], where it was shown that the map (|7.7(l is 
positive if and only if 



a>l, a + b + c>3, 1 < a < 2 =S> be > (2 - a) 2 



while it is decomposable if and only if 



a > 1, l<a<3 



be > (3 - a) 2 /4. 



(7.8) 



(7.9) 



In particular, ^2,o,c is positive but not decomposable for c > 1. All generalized indecomposable Choi maps are 
known to be atomic [33j |. it is they cannot be written as a convex sum of 2-positive and 2-co-positive maps |97|. 
An example of an indecomposable map belonging to V4 was given by Robertson 1861 . A family of indecomposable 
maps for an arbitrary finite dimension N > 3 was recently found by Kossakowski 57]. They consist of an affinc 
contraction of the set M.( N > of density matrices into the ball inscribed in it followed by a generic rotation from 
0{N 2 — 1). Although several other methods of constructio n of indeco mposable maps were proposed |5ll l78l l97l l98 | . 
some of them in the context of quantum entanglement llOOL IllOj , the general problem of describing all positive 
maps remains open. In particular, it is not known, if one could find a finite set of K positive maps {^j}, such that 

V N =conv hull (uf =1 ^j(CV N )) ■ 

Due to the theorem of St0rmer and Woronowicz the answer is known only for N = 2, for which K — 2, \&i — 1 and 
^2 = T. As emphasized by Horodeccy in an important work |37| these properties of the set Vn become decisive for 
the separability problem: the separability criterion based on the positivity of (1 ® T)p works for the system of two 
qubits, but does not solve the problem in the general case of the N x N composite system. 




decomposable maps ,** 

CP \^Ccf> 



non - decomposable maps 




a> V 2 = positive linear maps M m ^M. m b) T\ ^positive linear maps A^ (J) ^-Ai m 



FIG. 6: Sketch of the set positive maps: a) for N = 2 all maps are decomposable, b) for N > 2 there exist non-decomposable 
maps. 



The indecomposable maps are worth investigating, since each such map provides a criterion for separability. Con- 
ditions for a positive map $ to be decomposable were found some time ago by St0rmer 92]. Since this criterion is 
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not a constructive one, we describe here a simple test which may confirm the decomposability. Assume first that the 
map is not symmetric with respect to the transposition $ ^ T<£>. These two different points determine a line in the 
space of maps, parametrized by /3, along which we analyze the dynamical matrix 

£>/3*+(l-/3)T$ = + (1 - P)Dl A (7.10) 

and check its positivity by diagonalization. Assume that this matrix is found to be positive for some /?» < (or 
/3* > 1), then the line l|7.10|l crosses the set of completely positive maps (see Fig. E^). Since D(fl*) represents a 
CP map ^cp, hence D(l — /3*) defines a completely co-positive map 'fop, and we find an explicit decomposition, 
4? = [— P^cp + (1 — /3*) v I'ccp]/(1 — 2/3*). In this way the decomposability of 4? may be established, but with this 
criterion one cannot confirm that a given map is indecomposable. 




a) decomposable map <£> b) non -decomposable map <E> 



FIG. 7: Geometric criterion to verify decomposability of a map (a) if the line passing through $ and $T crosses the set of 
completely positive maps, a decomposition of <£> is explicitly constructed. 



To study the geometry of the set of positive maps it is convenient to work with the Hilbert-Schmidt distance, 
defined by the HS norm of the difference between the superoperators, d(^,$) = dHs{L^,L§) = \\L^ — L^\\hs- 
Since the reshuffling of a matrix does not influence its HS norm, the distance can be measured directly in the space 
of dynamical matrices, = cIhs{D^i,D^). Note that for unital one qubit maps, (|6.4() with k = 0, one has 

d(<E>i, <!>2) = 1 77*1 — f)2\, so Fig. OH represents correctly the HS geometry of the space of N = 2 unital maps. 

In order to characterize, to what extent a given map 4/ acting on MS N > is close to the boundary of the set of 
positive (CP or CcP) maps, let us introduce the following quantities 

a) complete positivity, cp($) = min (p\D$\p), (7-H) 

b) complete co — positivity, ccp($) = min {p\D^ A \p) , (7-12) 

c) positivity, p($) = min \(x <E> y\D$,\x (g) y)\. (7-13) 

x,y£CP N ~^ 

The two first quantities may be easily found by diagonalization, cp(<I>) = min{cig(Z)$)} and ccp(4>) = 
min{eig(D^ A )}. Although p(<&) > cp(4>) by construction the evaluation of positivity is more involved, since one 
needs to perform the minimization over the space of all product states, i.e. the Cartesian product CP^ -1 x CP^" 1 . 
No straightforward method of computing this minimum is kno wn, so one has to rely on numerical minimization. In 
certain cases this quantity was estimated analytically by Terhal |100| and numerically by Guhne et al. 0,H2 in the 
context of characterizing the entanglement witnesses. In fact a non-positive dynamical matrix D, which describes a 
non completely positive map, may be just considered as an entanglement witnesses - an oper ator D such TvDp is not 
negative for all separable states and negative for a given entangled state 0, |6^, |(5(| |83l Il0l| . 

As follows from the positivity of D and D Ta and the property of block positivity l|4.7|l . a given map $ is completely 
positive (CcP, positive) if and only if the complete positivity (ccp, positivity) is non-negative. As marked in Fig. [SJd, 
the relation cp($) = defines the boundary of the set CPn, while ccp(<I>) = and p(Q) = define the boundaries 
of CcP n and Vn- By direct diagonalization of the dynamical matrix we find that cp(l) = ccp(T) = and ccp(l) = 
cp(T) = -l. 

For any not completely positive map "^nCP one may look for its best approximation with a physically realizable CP 
map 4>cp, e.g. by minimizing their HS distance d(4?nCP, ^cp) _ see Fig. (5^. Such maps, called structural physical 
approximation were introduced in [4^ to propose an experimentally feasible scheme of entanglement detection and 
later studied in [28| . 

To see a simple application of complete positivity, consider a non physical positive map with cp($ n cp) = — x < 0. 
Its possible CP approximation, but generally not the optimal one, may be constructed out of its convex combination 
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TABLE V: Jamiolkowski Isomorphism 18.31 between Hermicity reserving linear maps $ definded on the space of mixed states 
which act on TLm and the operators D<s, on the composed space TLn ®Hn- 



Isomorphism 


Linear maps 


Hermitian operators 
-D$ : 7~L N 2 — > 7~t]\[2 


Ji 


positive maps $ 


operators D 
positive on product states 


Ju 


completely positive maps $ 


positive operators D 


r 

Jin 


quantum operations: 
CP, trace preserving maps 


states p = D/N such that 
TtaD = 1 


Jllla 


completely positive, unital 
maps 


states p = D/N such that 
Tr B D = 1 


Jiv 


CP and CcP 
trace preserving maps 


mixed states p = D/N 
such that TrA-D = 1 with 
positive partial transpose 


example 
of Jiv 
for N = 2 


CP and CcP 
trace preserving maps 
$ : M {2) -> M {2) 


separable states p — D/N 
of a 2 x 2 composite system 
such that TtaD = 1 


JlVa 


unitary rotations 
p' = UpU\ 
= {U® U*) R 


maximally entangled 
pure states 
(U® tW)(ip\ 


N = 2 example 

of Jiv a 
Pauli matrices 

versus 
Bell states 
P<t> = 10X01 


1^(1,0,0,1) 
a x <-» (0,1,1,0) 
<J» <-» (0, -i,i,0) 
a z ^ (1,0,0,-1) 


h/>+> = ^|00> + |11) 

|0 + ) = -^|Ol> + |lO> 
|0_> = ^|01> - |10) 
|^_> = ^|00>-|11> 


Jv 


completely depolarizing 
channel $» 


maximally mixed 
state p, = 1/N 



with the completely depolarizing channel 'J'*. Diagonalizing the dynamical matrix representing the map \t x = 
a3>nCP + (1 — a)^* with a = l/(Nx + 1) we see that its smallest eigenvalue is equal to zero, so fy x belongs to the 
boundary of CV n ■ Hence the distance <i($ n cp, 3> x )j which is a function of the complete positivity cp(<I> n cp), gives an 
upper bound for the distance of $ n cp from the set CP. In a similar way one may use ccP(<!>) to obtain an upper bound 
for the distance of an analyzed non CcP map $ n c c p from the set CcP . Interestingly, the solution of the analogous 
problem in the space of density matrices allows one to characterize the entanglement of a two-qubit mixed state p\ by 
its minimal distance to the set of separable states. In the two-qubit system the entangled states do not have positive 
partial transpose, so T^(pi) is not a state. As shown in Fig. Eh one may also study a related problem of finding the 
state p2 which is closest to Ta(pi). 



VIII. JAMIOLKOWSKI ISOMORPHISM 

Let CP n denote the space of all trace preserving, completely positive maps $ : MS N ^ — > A4^ N \ Note that this is a 
convex set. Any such a map <!> may be uniquely represented by its dynamical matrix D$ of size iV 2 . It is a positive, 
Hermitian matrix and its trace is equal to N. Hence the rescaled matrix /?$ = D$/N represents a mixed state in 
MS N \ and the entropy of the operation S{<&) equals to the von Neumann entropy of p$. In fact rescaled dynamical 
matrices explore only a subspace of this set determined by the trace preserving conditions (|4.6|) . which impose A^ 2 

constraints. Let us denote this A^ 4 — A^ 2 dimensional set by M.j . Since for any trace preserving CP map we may 

find a dynamical matrix, and vice versa, the correspondence between the maps from CP n and the states of M\ is 
one-to-one. In Table Ivl this isomorphism is labeled by J/77. 

Let us find the dynamical matrix for the identity operator, L = t N 2 , 

ijiu = ^mn^nv so that Dmfi = {Lmfi)^ = 8 m ^8 ni/ — N p m „, (8-1) 
where p^ = IV'XV'I represents the operator of projection on the maximally entangled state of the composite system 



N 1 
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This state is written in its Schmidt form (for TV = 2 it is the famous Bell state), and we see that all its Schmidt 
coefficients are equal, Ai = A.; = Ajy = 1/N. Thus we have found that the identity operator corresponds to the 
maximally entangled pure state \ip){^>\ of the composite system. Interestingly, this correspondence may be extended 
for other operations, or in general, for arbitrary linear maps. Any linear map $ acting on the space of mixed states 
M.^ N > can be associated, via its dynamical matrix with an operator acting in the enlarged Hilbert state Hn^Hn 



$ : M 



(N) 



M 



(AT) 



D$ = JV[$<8 i\(\i>)(ip\) 



(8.3) 



To show this we represent the linear map $ by its matrix L introduced in 1(4. 1() . write the operator $ (8) 1 as an 
eight-indices matrix and study its action on the state p^ expressed by ((8.1(1 . 



L 



mn 
m'n 



1 ^ Np v , , 

H'v' "*u 



Lmn — T)r t 
M" ? 



(8.4) 



An analogous operation 1 ® $ acting on p^ leads to the matrix D F with the same spectrum. Conversely, for 
any positive matrix D we find the corresponding map $ by diagonalization. The reshaped eigenvectors of D, 
rescaled by the roots of its eigenvalues give the canonical Kraus form l|4. 914.13(1 of the corresponding operation $. 
Furthermore, the entropy of the quantum operation S(<f>) equals the von Neumann entropy S(p$) of the corresponding 
state p$ = Dq,/N. 



|il>-t-><"i|; 




b) states 



FIG. 8: Duality between quantum map $ acting on a part of the maximally entangled state and the resulting density 

matrix p = jfDi> 



Consider now a more general case in which p denotes a state acting on a composite Hilbert space TLm ®"Hn ■ Let <J> 
be an arbitrary map which sends into itself and let D%, = L R denote its dynamical matrix (of size TV 2 ). Acting 

with the extended map on p we find its image p' = [<£> £g> l](/o). Writing down the explicit form of the corresponding 
linear map in analogy to 1(8. 4J1 . and contracting over the four indices representing 1, we obtain 

(p') R = L< iP R so that p 1 = (D R p R ) R . (8.5) 

In the above formula the standard multiplication of square matrices takes place, in contrast to Eq. 1(4.1(1 in which the 
state p acts on a simple Hilbert space and is treated as a vector. 

Note that Eq. may be obtained as a special case of ((8.5(1 if one takes for p the maximally entangled state 

((8.2|) . for which (p^) R — 1. Formula ((8.5(1 provides a useful application of the dynamical matrix corresponding to a 
map $ acting on a subsystem. Since the normalization of matrices does not influence positivity, this result implies 
the reshuffling lemma (17.511 . 

Formula 1(8. 3|) may also be used to find operators D associated with positive maps $ which are neither trace 
preserving nor complete positive. The correspondence between the set of positive linear maps and dynamical matrices 
acting in the composite space and positive on product states is called Jamiolkowski isomorphism since it follows 
from his results obtained in |4^. Some aspects of the duality between maps and states were recently investigated in 
0, |3(| . Let us mention here explicitly certain special cases of this isomorphism labeled by Ji in Table [V] The set 
of all completely positive maps $ is isomorphic to the set of all positive matrices D, (case Ju). Unital CP maps 
are related to dynamical matrices which satisfy trsD = 1 (case Jin a )- The set of all quantum operation, (i.e. the 
trace preserving, CP maps) corresponds to the set of positive matrices D fulfilling another constraint, tr^Z) = 1 - 
see Table [3 item Jm. An apparent asymmetry between the role of both subsystems is due to the particular choice 
of the relation 1(8.3(1 : if the operator 1 (g) $ is used instead, the subsystems A and B in the partial trace constraints 
need to be interchanged. 

An important case J\y of the isomorphism concerns the states with positive partial transpose, p TjL > 0. called 
briefly PPT states. Another case, J/y a , relates the set of unitary rotations, p' = < E'(p) = UpU* with the maximally 
entangled states, {U <g> l)|i/')(V ; l- The local unitary operation (U <g> 1) preserves the purity of a state and its Schmidt 
coefficients. Thus the set of unitary matrices U of size N is isomorphic to the set of the maximally entangled pure 
states of the composite N x N system. In particular, vectors obtained by reshaping the Pauli matrices ex; represent 
the Bell states in the computational basis, as listed in Table IVl Eventually, case Jy consists of a single, distinguished 
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point in both spaces: the completely depolarizing channel and the corresponding maximally mixed state p». Note 
the following inclusion relations of the sets mentioned in Table [V] Jj D Jjj D Jm D Jjv 3 Jy C Jin a and 
J II D J nia D J iv a C J///, as sketched in Fig. El 



IX. QUANTUM MAPS AND QUANTUM STATES 



Relation (|8.3[) allows one to link an arbitrary linear map $ with the corresponding linear operators given by the 
dynamical matrix D$. Expressing the maximally entangled state in (|8.3I) by its Schmidt form l|8.2[l we may 
compute the matrix elements of D§ in the product basis consisting of the states \i <g> j). Due to the factorization of 
the right hand side we see that the double sum describing p$ = D$ /N drops out and the result reads 



(fc®»iD»ij<8i) = (k\mi)(j\)\i). 



(9.1) 



This equation may also be understood as a definition of a map $ related to the linear operator It proves the 
isomorphism Ji from Table if Z)$ is block positive, then the corresponding map $ sends positive operators into 
positive operators 48] . 

As listed in Table M and shown in Fig. ED the Jamiolkowski isomorphism (|9.1|) may be applied in various setups 0. 
Relating linear maps from Vn with the operators acting on an extended space TLm ® 7~Ln we may compare: 

i) individual objects, e.g. completely depolarizing channel and the maximally mixed state p*, 

ii) families of objects, e.g. the depolarizing channels and the Werner states |lf)5j . 

iii) entire sets, e.g. the set of CV U CcP maps and the set of PPT states, 

iv) certain problems, e.g. for an arbitrary CP map $ find the closest CcP map, and the problem of finding the 
PPT state (separable state for N = 2) closest to an arbitrary (entangled) state, and 

v) their solutions... 

To get some more experience concerning the analyzed duality between quantum maps and quantum states compare 
both sides of Fig. [5J Note that this illustration may also be considered as a strict representation of a fragment of the 
space of one-qubit unital maps (a) or the space of two-qubits density matrices in the HS geometry (b) . It is nothing 
else but the cross-section of the cube representing the positive maps in Fig. |3t along the plane determined by 1,T 
and $*. 

T A W + > <^ + |) 

"i 

family of 
depolarising 
maps 




^positive linear maps M 



b) Hermitian operators acting in H. 



FIG. 9: Isomorphism between objects, sets, and problems: a) linear one qubit maps, b) linear operators acting in two-qubit 
Hilbert space TI4. Labels Ji refer to the sets related by the isomorphisms defined in Table [V] 



Let us finish this section by pointing out an analogue of the Jamiolkowski isomorphism in the classical case. The 
space of all classical states - probability vectors of size N - forms the N — 1 dimensional simplex Ajv-i- A discrete 
dynamics in this space is given by a stochastic transition matrix Tjy : Ajy-i — * Ajy-i- It contains non-negative 
entries, and due to stochasticity the sum of all its elements is equal to N. Hence reshaping the transition matrix 
Tjv and rescaling it by 1/N we receive a probability vector t of length N 2 . The classical states defined in this way 
form a measure zero, N(N — 1) dimensional, convex subset of A N 2_ 1 . Consider, for instance, the set of N = 2 

a b 
1-a 1-6 



stochastic matrices, which can be parametrized as T2 



with a, b € [0, 1]. The set of the corresponding 
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probability vectors t = (a, 6, 1 — a, 1 — 6)/2 forms a square of size 1/2 - the maximal square which may be inscribed 
into the unit tetrahedron A 3 of all N — 4 probability vectors. 

The classical dynamics may be considered as a (very) special subclass of quantum dynamics, defined on the set 
of diagonal density matrices. Hence the classical and quantum duality between maps and states may be succinctly 
summarized in the following, commutative diagram 

quantum: [$ : M {N) -> M {N) ] — > G M {n2) 

i maps | states (9.2) 

classical : [T : Ajv-i — > Ajv-i] — ► t* G Ajv2_i . 

Alternatively, vertical arrows may be interpreted as the action of the coarse graining operation \&cg defined in section 
IVI For instance, for the trivial (do nothing) one-qubit quantum map the superopcrator L = I4 restricted to 
diagonal matrices gives the identity matrix, T = 1 2 , and the classical state t = (1, 0, 0, l)/2 G A 3 . On the other hand, 
this very vector represents the diagonal of the maximally entangled state \D$, = \ip)(ip\. To prove the commutativity 
of the diagram (|9.2() in the general case define the stochastic matrix T as a submatrix of the superoperator Q4.1[l . 
T m n = Lmm (left vertical arrow). Note that the vector t obtained by its reshaping satisfies t = diag(L ii ) = diag(-D$), 
so it represents the diagonal of the dynamical matrix (right arrow). 



X. ENVOI 



In this work we analyzed the space Vn of positive maps, which send the set of mixed quantum states into 
itself, and its various subsets: the set CVn of completely positive maps, and the set Bn of bistochastic maps. For 
any quantum map \& we have introduced three quantities, H7.lll7.13fl . to characterize the location of '5 with respect 
to the boundaries of the sets of positive (completely positive or completely co-positive) maps. We have defined and 
investigated the set of unistochastic maps, the elements of which are determined by a unitary matrix U of size N 2 
and correspond to the coupling with the TV-dimensional environment, initially in the maximally mixed state. 

The spaces of quantum maps (stochastic, bistochastic, unitary evolutions), which act on the set A4^ N ^ of mixed 
quantum states, may be related with the corresponding classical maps (stochastic, bistochastic, permutations), which 
act on the simplex An_i of all TV-point probability distributions. We have extended the Jamiolkowski isomorphism, 
linking the space of quantum maps with the space of bipartite quantum states, to the classical case (|9.2|l . 

In general the state-map duality, analyzed in this paper, may be formulated and applied in a variety of contexts 
and set ups. For instance, the set of SU(i) matrices may be considered as: 

a) the space of maximally entangled states of a composite, 4x4 syste m, G CP 15 C A4^ 16 ' , 

b) the set of two-qubit unitary gates acting on - see e.g. |3a l58l l75l Ill3| . 

c) the set BU2 of one qubit unistochastic operations (|3.10(l . G BU2 C 62 C CV2- 

To conclude, we would like to convey two main points. On one hand, the set of positive maps has an interesting 
geometry, which is worthy of investigation in the general case N > 3. On the other hand, the marvelous duality 
between quantum maps and quantum states, based on the Jamiolkowski isomorphism, allows one to apply all the 
knowledge gained in studying the set of quantum maps to describe the set of quantum states (or conversely). For 
instance, the known structure of the set of completely positive maps relates to the structure of all quantum states, 
while the description of the (larger) set of all positive maps would allow us to describe the subset of all separable 
states. 
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M. Kus, W. A. Majewski and F. Mintert for fruitful discussions and helpful comments. Financial support by Komitet 
Badah Naukowych under the grant PBZ-Min-008/P03/03 and Swedish grant from VR is gratefully acknowledged. 
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